Condensate superfluidity and infrared structure: the Josephson relation 
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Although the phenomena of superfluidity and Rose 
condensation [1] are intimately related, the connection 
between the condensate density, no, and the superfluid 
mass density, ps, is only indirect. In the limit of di- 
lute gases the superfluid mass density simply equals the 
condensate density times the particle mass, m. More 
generally, the two quantities are related via the infrared 
structure of the single particle Green's function, G. The 
detailed relation was given by Josephson 0, [3| , 



(1) 



where G'ii(fc,0) is the single particle Green's function at 
momentum k and zero Matsubara frequency. 

This relation, applied in the critical region of Bose con- 
densation in a three dimensional system, determines the 
scaling structure of the superfluid density Further- 
more, it implies that in a two dimensional interacting 
Bose gas, superfluidity can exist below the Kosterlitz- 
Thouless-Berezinskii transition temperature, even in the 
absence of condensation in the thermodynamic limit Q . 
The Josephson relation was derived originally by heuris- 
tic arguments. In this note we derive the Josephson re- 
lation at finite temperature within the framework of di- 
agrammatic perturbation theory, for infinite as well as 
large but finite size systems. The present derivation, 
which follows a different path from that of Josephson, 
extends to finite temperatures the discussion of Gavoret 
and Nozieres Q of the microscopic connection between 
superfiuidity and condensation at zero temperature. 

A superfluid effectively has an extra dynamical degree 
of freedom, the velocity of the superfluid, Vs, with respect 
to the wall velocity [3]. In the presence of superflow the 
free energy density, F, of the system is 



Fivs,T,p) 



F{0,T,p) + ^p,vl 



(2) 



where T is the temperature and p the chemical potential. 
Thus thermodynamically. 



Ps 



F{vs,T, p)\v,=o- 



(3) 



and the normal density p„ is given by pn = p — Ps where 
p is the mass density of the system. 



The Hamiltonian H{vs) in the rest frame of the super- 
fluid, written in terms of the variables of the "lab" frame, 
in which the walls are at rest, is: 



H{vs) =H- P,vs 
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(4) 



where P is the total momentum operator. We work in a 
system of large volume V , and assume periodic bound- 
ary conditions in the z direction fs']. We consider the 
free energy density of a system with non-zero superfluid 
velocity with respect to the walls, in the frame moving 
with the superfluid. Differentiating the partition func- 
tion with respect to Vs^ we obtain 



dFldv, = ^{P,-Mv,)/V, 



(5) 



and then, noting that the total momentum operator in 
the z direction commutes with H, we have 



d^F/dvl = p~l3{PjMV, 



(6) 



'At Vs — Q. Thus the normal mass density is given by 

P„=/3(P2)/K (7) 

In a normal system, the total momentum P has a 
classical Boltzmann distribution with probability oc 
exp(— /3P^/2M), where M is the total mass, which im- 
plies I3{P'^)/V = M/V = p, and thus /?„ = p. Only devi- 
ations of the distribution of states with total momentum 
P from classical can give rise to superfluidity; in the su- 
perfluid phase, the distribution of total momentum is no 
longer classical as a consequence of entanglement of the 
total momentum and superfluid velocity. 

Microscopically the superfluid and normal mass densi- 
ties are defined in terms of the transverse current-current 
correlation function 0, Q. In fact, the microscopic def- 
inition, is, as we shall see, equivalent to Eq. d?]). We 
consider a superfluid contained in an inflnitely long cylin- 
drical vessel oriented along the z axis. If the walls of the 
container move with a small velocity in the z direction, 
the normal mass density moves with the walls, while the 
superfluid remains at rest in the laboratory frame. The 
response of the fluid flow to motion of the walls is speci- 
fled by the current-current correlation function 

T„-(r,r',L.)= /die^-(*-*')([j,(r,t),j,(r',0]>, (8) 
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where the current density j is given by 

j(r) = ^ [V-t^V^W - V(r)] , (9) 

with tpir) the particle annihilation operator. The Fourier 
transform to momentum k and real frequency u can 
be decomposed into longitudinal (L) and transverse (T) 
parts as 



fc2 



The longitudinal and transverse components of the time- 
ordered current-current correlation function in imaginary 
time have Fourier transforms to Matsubara frequency , 



27r — LO 



ill) 



The long-wavelength static longitudinal response func- 
tion - describing the response of the system in a tube 
with closed ends moved along the z direction [1, Q - 
obeys 

SL(fc^O,0)= lim lim Sz2(k, = 0) = -p/m^, 

— ^0 — >0 

(12) 

where _L denotes the directions perpendicular to z. Then 

duj T ^(k, uj) 



p = lim 

fc^O 



27r 



(13) 



where p is the total mass density. This relation follows 
from the f-sum rule together with the continuity equa- 
tion. On the other hand, the long wavelength limit of the 
transverse response, which describes the response of the 
system in an open ended moving container (or with pe- 
riodic boundary conditions), reduces to the normal mass 
density, p^ = p - p^: 



^T{k ^ 0,0) = lim lim ^2^(k, z,, = 0) = -pn/m 

k I ^0 kz^Q 



lim TO^ 

fe^O 



duj Tx^k, u) 
27r UJ 



(14) 



(15) 



While the order of limits in Eq. p2)l describes a system 
closed in the z direction, the order in Eq. (fT4|) describes 
a system open in the z direction [3i]. 

A superfluid is characterized by pn < p- As long as the 
two-particle Green's function is regular, the kx,ky,kz — *■ 
limits in Eqs. (fT2|) and (fT4|) cannot depend on the order 
of limits, and thus p„ — p. Superfluidity implies that the 
two-particle Green's function behaves singularly in the 
infrared limit. 

Furthermore, since / drj{rt) = P{t)/m, we have 



lim lim 0) 



-03 



dt{p,{t)PM), 



(16) 



where V is the system volume. Since for an infinite sys- 
tem, Pz is independent of we rederive Eq. ([7]). 

Having established the equivalence of the two defini- 
tions of ps , we turn to deriving Josephson's relation in a 
Bose system, using the definition ([3]). We take the free 
energy density to be a functional of the order parame- 
ter, (iplr)), of the system and the single particle matrix 
Green's function, 

Girt,r't') = 
-i((r(*(rt)*t(i.Y))) - (*t(i.'t'))(*(rt))) . (17) 

Here the two-component field operator is ^'(rt) — 
(i/'(rt), ^^(rt)). In equilibrium, the first order varia- 
tion of F with respect to the order parameter vanishes, 
(5F/(5(*) = SF/Si'^/)* = 0. Static variations of F with 
respect to the order parameter are then given in terms of 
the inverse of the single particle matrix Green's function, 

VSF^-^ y"drdr'5(*(r))*g-i(r,r')(5(*(r')), (18) 

where (^P) = (V')*) is the two-component order pa- 
rameter. 

The dependence of F on the superfiuid velocity enters 
through the phase of the order parameter. 



(19) 



where rip is the condensate density evaluated at = 0. 
To second order in Vs , the dependence of the free energy 
density on the superfiuid velocity arises solely from the 
momentum mvs of the condensate wave function. To 
lowest order in Vs, S{'ip{r)) = imVs ■ r^/no, so that the 
second order change in F is, 



SF = 



^ y"drdr'(v, •r)(v, -r') 

y<{ig-')ii-ig-')i2) (r,r') 



lim (v, • Vi.) 

X ((^-^)ii(fc,0)-(g-i)i2(fc,0)) 



(20) 

Since the inverse of the static Green's function has the 
form, 



— 2j21 p. — £k~ 2^22 



(21) 



where = k? j^m and the (fc, 0) are the self-energies, 
the variation of the free energy density (PO)) is given in 
terms of the S's by 

1 c)'^ 

5F = -nom2i;2_(£^ + l],,(fc,0)-Si2(fc,0))^,^Q. 

(22) 



3 



After brief algebra we obtain, 

(Gii)-^ (fc,0) = -2 [efe - A* + Sn(fc,0) - Ei2(fc,0)] 

[A^-efc^S]n(fc,0) + Si2(fc,0)]' 
/i- Efc - Eii(A;,0) 

(23) 

We assume that below the critical temperature, Tc, 
the single particle excitation spectrum is gapless. Then 
Eq. reduces to Eq. ((T|). In this case, the chemical 
potential, /i, which depends on tiq, is determined by the 
Hugenholtz-Pines relation 

A* = £11(0,0)- £12(0,0). (24) 

Thus the final term in Eq. (|23p is of order k'^, and can 
be neglected as fc ^ 0. Rewriting Eq. ((22|) in terms of 
G^^{k,0), we find 

1 9^ 1 

which implies the Josephson relation ([1]). 

Equivalently, we can derive p„ by carrying out a gauge 
transformation on the G^- in which 

Gii(r,r',0) ^ e™^" ("^-'-')Gii(r,r',0), 

Gi2(r,r',0) e™^"-("-+"^')Gi2(r,r',0), 

G2i(r,r',0) ^ e-™^"-(--+-^')G2i(r,r',0), 

G22(r,r',0) ^ e-™^"-("-'-')G22(r,r',0). (26) 

Under this transformation, the free energy density 
changes by 

SF ^ ^p,,vl. (27) 

Now the first order variation of the free energy density 
with SG gives a purely kinetic contribution, 

<5Ffc,„ - ^l.Tr{g-'Sg) (28) 

= ^(/?-mrio)w^, (29) 

since (^o-i)ii(r,r') = (a* + V2/2m)^(r - r'), 
(e„-i)ii(r,r') ^ (go-^)22(r,r'), and (go"i)i2(r, r') = 



(5g^_[2i(r, r') = 0. The second order variation is given 

by M 

SF„u = - 2^ TriTr2(5g(l)£-i(l, 2)5^(2). (30) 

The connected two particle correlation function, 
£(1, 2) = ^^2(1, 2) - ^^(1)^^(2), is the two particle Green's 
function, t/2, with the uncorrelated one particle Green's 
function contribution subtracted. To evaluate Eq. ((30|) 
with we note that the under the gauge trans- 

formation coupled with the transformation {ip{r)) — > 
(■0(1")) e™''^" '', the diagrams involving the interaction are 
invariant. Thus from Eq. (P^ . we see that ([50)1 con- 
tributes 

d 

SFint = -no^Ti^w^^p (Sii(fc,0) - Ei2(fc,0))j,^o ■ 

(31) 

Assembling the pieces, and using Eq. ((25)) we find that 
under the transformation (|26)) . 

SF = ^(^p~ nom' lim -^G^,\k, 0)) 

= ^(P-P.))^', (32) 
as expected. 

We note that our derivation of the Josephson relation 
is in fact valid in a finite size system. There the Green's 
functions are defined only on a discrete set of points in k 
space, the limit A; ^ in Eq. ([T]) is replaced by the limit 
k ^ kfj with /co = 27r/L, and the Josephson relation, be- 
comes discretized. In addition, the relation remains valid 
in two dimensions in the gapless phase 0], since we did 
not need to explicitly specify the number of dimensions 
in the derivation. 
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